Reduction of velocity fluctuations in a turbulent flow of gallium by an external 

magnetic fleld. 
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The magnetic field of planets or stars is generated by the motion of a conducting fluid through a 
dynamo instability. The saturation of the magnetic field occurs through the reaction of the Lorentz 
force on the flow. In relation to this phenomenon, we study the effect of a magnetic fleld on a 
turbulent flow of liquid Gallium. The measurement of electric potential differences provides a signal 
related to the local velocity fluctuations. We observe a reduction of velocity fluctuations at all 
frequencies in the spectrum when the magnetic fleld is increased. 
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Magnetic fields of planets and stars are created by the 
dynamo instability, which converts part of the kinetic en- 
ergy of an electrically conducting fluid into electromag- 
netic energy. Recently, several groups succeeded in gen- 
erating fluid dynamos using liquid sodium flows [U, U Q . 
The feedback of the magnetic field on the flow through 
the Lorentz force causes the saturation of this instabil- 
ity. This effect is quite well understood in the Riga Q 
and Karlsruhe [2] dynamos in terms of mean field mag- 
netohydrodynamics or of the change of the velocity pro- 
files 0, I3j [1, @|- In these two cases the turbulence in- 
tensity is low and the turbulent motions are restricted to 
small scales. In the case of the VKS dynamo [l| for exam- 
ple, the turbulent fluctuations are present at all scales of 
the experiment. We are interested in the detailed dynam- 
ical mechanism involved in this feedback on the turbulent 
flow, which remains poorly understood. In particular, 
it is not clear whether the magnetic field first inhibits 
the smallest turbulent eddies, i.e. the ones containing 
the lowest amount of kinetic energy, as first assumed 
by Batchelor Q, or if it affects all scales of the flow. 
Several experiments studied the influence of strong mag- 
netic fields on decaying turbulent fiows of liquid metal 
d, 1^, [13, [iH ■ Most of them concern freely decaying tur- 
bulence and are performed in a range of dimensionless 
parameters very far from the one involved for weakly su- 
percritical dynamos. Here we focus on the braking effect 
of the magnetic fleld on a stationary turbulent flow. A 
constant magnetic field of tunable amplitude is applied to 
a gallium fiow. The statistics of the velocity field can be 
accessed through the measurement of electrical potential 
differences in the metal flow. We report the decrease in 
the intensity of the velocity fluctuations at all frequencies 
as the magnetic field is increased. 

One liter of gallium fills a glass cylindrical cell, 12.5 cm 
in diameter and 7 cm in height (fig[l|). A coil imposes a 
constant vertical magnetic field that can be varied up to 
1 kG. A propeller is positioned 3 cm off- axis and 3 cm 
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FIG. 1: Cross-section of the experimental setup. A propeller 
is placed off axis in a cylindrical cell flUed with liquid gallium. 
The vertical magnetic fleld is generated by a cylindrical coil 
and is monitored with a gaussmeter. An electrical potential 
difference is measured by a pair of copper electrodes. 



above the bottom of the vessel. Its radius is i? = 3.5 cm. 
The propeller is rotated so that it ejects the fluid towards 
the bottom of the cell. The off centering of the propeller 
prevents the fluid from rotating as a solid body. This 
generates a highly fluctuating flow. Rotation frequencies 
frot of the propeller range from 2 to 7 Hz. 

In the approximations of magnctohydrodynamics, the 
magnetic field B follows the induction equation: 



OB 

— = V X (v X B) + ?7AB , 



(1) 



where rj = {figa)^^ is the magnetic diffusivity {a is the 
electrical conductivity) and v is the flow velocity. For an 
incompressible fluid, the latter is governed by the Navier- 
Stokes equation with the Lorentz force: 
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V)v = - Vp + pi^Av + j X B , 



(2) 
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where v is the kinematic viscosity, p is the pressure, j is 
the electrical current density and p is the density of the 
fluid. Three dimensionless numbers can be defined: (i) 
the usual hydrodynamic Reynolds number: Re = LV/v 
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(V and L being velocity and length scales) , (ii) the mag- 
netic Reynolds number = LV/rj which is the ratio 
of the induction term over the magnetic diffusive term 
(Ohmic diffusion) in the induction equation, (iii) the in- 
teraction parameter N which compares the Lorcntz force 
j X B to the convective acceleration in the Navier-Stokes 
equation. For an applied magnetic field Bq, the induced 
current is j '--^ (t(v x Bq), so that N = aLBo^ / pV 

In the case of gallium one has a = 3.8610^ il^^m"^, 

^ 3.11 10"'^ m^s-i and p = 6090 kg/m^. 

In our experiment, the Reynolds number can be var- 
ied in the range SIO'' < Re < 1.7 10^ so that the flow 
is turbulent. The magnetic Reynolds number is in the 
range 0.07 < Rm < 0.25. The interaction parameter can 
be rewritten N — aBo^ /{2nfrotp) by using R frot a 
velocity scale. 

Velocimetry in liquid metals can be achieved by mea- 
suring the electromotive force that arises from the motion 
of the metal in an applied magnetic field. Measurements 
of the voltage between two electrodes, a few millimeters 
apart, have been performed for more than 50 years. The 
magnetic field can either be a large scale field [111. [l^ . [isj 
or the localized field of a small magnet [3, [lB| . We use 
the first configuration in our experiment. The electrodes 
are located at the end of a cylindrical rod plunging into 
the liquid metal (fig. [T|). Each one is made of a 5 mm 
long copper wire, insulated except at the very end which 
is in electrical contact with the fluid. They are separated 
by a 3.5 mm wide gap. 

In the approximations of magnetohydrodynamics. 
Ohm's law in a moving fluid is written : j = (^(E+vxB). 
In quasi-static regime one can consider that V • j — 0. 
With the choice V • A = 0, the divergence of Ohm's law 
yields : 

A(/) = • B - ^ov • j (3) 

The potentials (f) and A are defined by E = — V0 — ^ 
and Lo is the vorticity of the flow. The order of magnitude 
of the right-hand side terms of ((3|) are : cu ■ Bq ~ BqV/Ik 
and /iqv • j ~ hV /la- b is the order of magnitude of 
the magnetic field induced from Bq by the fiuid mo- 
tion. Ik ~ LRe~^^^ is the Kolmogorov length (dissi- 
pative length scale of the velocity) and ^ LRmT^^'^ is 
the length scale of ohmic dissipation of the induced cur- 
rents. L is chosen as the integral length scale of the flow. 
In our experiment, Rm is low so b ^ Bq and the second 
r.h.s term is negligible so that 

A(I) = lu-Bq (4) 

Thus potential measurements provide a measurement of 
the component of vorticity parallel to the applied held. 
Similarly to the case of pressure, a direct measurement 
of (j) would be nonlocal, because of the Laplacian oper- 
ator. For length scales larger than the separation I of 
the electrodes, the measurement of potential differences 
gives the gradient of 0. A signal of order vBq is ob- 
tained (one spatial integration of ((3|)). Thus, we expect 



the potential difference measurements to be related to 
the dynamics of the velocity fluctuations, but a precise 
calibration of the probe is not straightforward. For the 
fluctuations corresponding to length scales much smaller 
than /, the difference between two decorrelated values of 
the potential (j) is measured. We expect the statistics of 
the potential difference to be the same as that of the po- 
tential itself. In the framework of Kolmogorov's theory of 
turbulence, the Fourier spectrum of vorticity is supposed 
to behave as k^^^. Because of the sweeping of velocity 
structures by the turbulent flow or the Taylor hypothesis 
(when relevant [3]), one gets directly a prediction for 
the frequency spectrum. We expect the spectrum of the 
potential difference to behave as f~^^^ for frequencies 
lower than a cutoff frequency fc V/2'kI and as /"^^^^ 
for higher frequencies. 
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FIG. 2: Top: temporal signal measured with the potential 
probe for a rotation frequency of propeller of 7 Hz under a 
magnetic field of 738 G. Bottom: corresponding power spec- 
trum. 

For a magnetic field of about 50 Gauss, a signal of a 
few microvolts is measured. It increases with both the 
applied field and the rotation frequency of the propeller. 
A time series of the potential is displayed in fig. [Altogether 
with the corresponding temporal power spectrum. The 
qualitative behavior of the spectrum is in agreement with 
the discussion of the previous paragraph. For intermedi- 
ate values of the frequencies, starting from the rotation 
frequency of the propeller, the spectrum decays following 
roughly a power law which exponent is close to —5/3. For 
frequencies larger than 40 Hz, a steeper power law decay 
is observed, with an exponent close to —11/3. Beyond 
100 Hz, the signal falls under the noise level. This is 
consistent with a hydrodynamical interpretation of our 
potential measurements. For lower driving speeds, the 
/~^/^ zone is not very visible, most likely because of the 
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FIG. 3: Standard deviation of the fluctuations of tfie potential 
difference divided by the rotation frequency of the propeller 
versus applied magnetic field . 
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FIG. 4: Standard deviation of dimensionless velocity fluctua- 
tions, versus interaction parameter N. 



lower value of the Reynolds number that may not de- 
velop a clear inertial range. The quantitative behavior of 
the spectra depends somewhat on the position and the 
orientation of the probe, because of the inhomogeneity of 
the flow. 

We evaluate quantitatively the amplitude of the fluc- 
tuations of the potential difference Scj) by the standard 
deviation cr^^ of the recorded signal. As seen in fig. [31 
when the applied magnetic field is gradually increased, 
the amplitude of 6(j) first grows linearly with Bq and frot- 
For larger values of the magnetic field, a departure from 
the linear behavior in Bq is observed which we ascribe to 
the effect of the Lorentz force on the fiow. 

The potential difference is divided by Bq (constant dur- 
ing a given measurement) and also by the width I of the 
probe. This quantity v — 5(j)/ B^l is homogeneous to a 
velocity and is related to the velocity fluctuations in the 




FIG. 5: Comparison between power spectra of velocity for 
a rotation frequency of 5 Hz for two different values of the 
external magnetic field. Inset : ratio between the two power 
spectrum densities. 



vicinity of the probe. Although the link is not direct be- 
tween velocity and potential difference, for simplicity, we 
call V (somewhat abusively) "velocity" in the following. 
The standard deviation of this velocity is non dimcn- 
sionalized by Rfrot and is noted <7y<- — (Ty/Rfrot in the 
following. Within the accuracy of our measurements, it 
decreases linearly with the interaction parameter, as seen 
in fig. m The decay of a^* can reach a factor 2. 

The obtained curves look similar to the potential mea- 
surements at large scale performed by Steenbeck et al. 
[13] ■ They recorded the potential produced by a effect 
at the border of a dedicated device with sodium flow. 
In their experiment, the average potential decreases with 
the applied field, and follows roughly a 1/N decay law. 
Their measurements are interpreted as a quenching of 
the a effect, corresponding to a braking of the flow. This 
saturates the efficiency of the induction processes and is 
responsible for the saturation of the a-dynamo ^ISj . It 
supports the model of the a effect, where the average 
electromotive force is written as (v x B) = aBo. Here 
we observe directly the reduction of the turbulent veloc- 
ity fiuctuations resulting from the magnetic braking of 
the fiow. 

One open question on the spectral behavior of the mag- 
netic braking phenomenon is whether it acts on a privi- 
leged scale. In particular, one could expect the braking 
to be more efficient at the smallest scales of the flow 
and therefore to damp preferentially the highest frequen- 
cies of the spectrum. Indeed, following Batchelor Q, a 
scale S would be damped if the volumic kinetic energy 
^pvg at this scale is below the volumic magnetic energy 
^Bo^/fiQ. The braking effect of the magnetic held could 
then be modeled as an enhanced dissipation. To answer 
this question, we compare the spectra of dcfi/Bo for two 
values of the applied magnetic field at a given frequency 
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frot of the propeller (fig. [5]). The spectrum corresponding 
to the highest Bq is below the other one but no relative 
variation in the damping is observed over the frequen- 
cies. The two spectra can be obtained one from another 
by a simple scaling factor which is the ratio of the two 
squared standard deviations cr^. . The magnetic braking 
is performed similarly at all scales in the present range 
of dimensionless parameters. 

Several authors [1, [1, [l^, [ll| studied the influence of a 
magnetic field on a flow both experimentally and numer- 
ically. For decaying turbulence, they observed a faster 
damping of the velocity field for large values of the inter- 
action parameter N. This occurs simultaneously with a 
bidimensionnalisation of the flow. In our experiment, no 
evolution is observed toward a f~'^ decay of the spectrum, 
expected for 2D turbulence. One difference may be due to 
the continuous forcing in our experiment but also to the 
different values of N and Rm. In this case the large scale 
eddy turnover time may be lower than the decay time of 
the velocity field via the Lorentz force, for moderate val- 
ues of N. The mechanical energy injection eo ^ pV^/L 
should be compared to the power of the Lorentz force. 
The latter is identified with the dissipation by Joule ef- 
fect j • E aV^Bo'^. The ratio of the two is N which 
takes values lower than one in our case. The dissipation 
of the kinetic energy is thus done primarily via the tur- 
bulent cascade at TV < 1. The current density is of order 
j ^ (tV X Bo, so that the Lorentz force is aBo^V. Thus 
we can build a braking time — p/aBo^, which must 
be compared to the eddy turnover time t/j - after 
which the velocity field is decorrelated. Once again, for N 
smaller than unity, one has Th/rm = N < 1, which means 
that the magnetic field does not have enough time to af- 
fect the structure of the fiow and to make it anisotropic. 
Consequently in the range of parameters of our experi- 
ment and also for the saturation of the dynamo effect. 



the application of a magnetic field to a turbulent flow, 
brings an additional braking term, without changing the 
spectral behavior of velocity fluctuations. On the con- 
trary, for TV > 1, a variation of the slope of the velocity 
spectrum with th e ap plied field is expected - as observed 
by Eckert et al. The total dissipation of the fiow 

increases with the magnetic field because of Joule dis- 
sipation. A smaller part of the injected power is avail- 
able for the turbulent cascade. Nevertheless it neither 
changes much of its dynamics nor breaks its scale invari- 
ance. An efficiency factor F < 1 depending on the inter- 
action parameter can be introduced to estimate the effec- 
tive kinetic dissipation rate. From the usual arguments 
of Kolmogorov's 1941 theory, the velocity spectrum can 

then be written E{k) ~ {V^L-^F{N)f'^ k-^/^ . Conse- 
quently when the magnetic field is increased, the velocity 
spectrum preserves its shape, but its amplitude decreases 
as F{Nf/^. 

In summary, we observed the influence of the Lorentz 
force on a turbulent fiow of gallium through measure- 
ments of electrical potential differences. A strong re- 
duction of the turbulent fluctuations is evidenced as the 
applied magnetic field is increased. In the range of in- 
teraction parameter and magnetic Reynolds number cor- 
responding to our measurements, the observed spectra 
show that the dynamics of turbulence is affected by the 
presence of the magnetic field only through a scaling ef- 
ficiency factor. 
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